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Abstract 

O" 

' Results for recoil nucleon induced polarization for (e, e'p) are presented using 

PsJ . various approximations for the relativistic nucleonic current, at the kinematics 

of a recent experiment at Bates. We see that the dynamical relativistic effects 
improve the agreement with the data. We make predictions for the induced 
Q> ' normal polarization and responses for TJNAF 89-033 and Mainz Al/2-93 

experiments. 
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!>. ■ I. INTRODUCTION 

o ■ 

Q , Several experiments have been proposed or have been carried out to measure the polar- 

ization of the ejected nucleon in (e, e'p) reactions In this way, new sets of polarization 

response functions can be isolated 0-0] • 
^ ■ If we write the cross section for the coincidence (e, e'p) reaction in terms of recoil nucleon 

2 ■ polarization dependent and independent terms, we have 



1 + P-a 
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where is the scattered electron energy, ctq is the unpolarized cross section, s denotes 
the nucleon spin projection upon a, and P is the induced polarization. Each of these 
observables can be written in terms of response functions that are bilinear combinations of 
the nuclear electromagnetic current operator [^^,0. If the electron beam is unpolarized and 
the experiment is performed in coplanar kinematics {(j)' = 0,7r), the relationship between 
the nuclear responses and the cross section is given by: 



rfV, E'\P' 
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Mott 

[Vl{Rl + RlSn) + Vt{Rt + RTSn) + 

COS <j)'VTL{RTL + RtlSu) + COS 20VTT(i?TT + RttSu)} 



The kinematical factors are Vl = Vr = A/2 + tan^ 6*6/2, Vtt = A/2, Vtl 



\J\ + tan^ 9e/2 and A = 1 — {uj/\q\)'^ where u and q are the energy and momentum transfer 
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in the reaction, 9^. is the electron scattering angle and E' , \P'\ are the energy and momentum 
of the ejected nucleon. Hence, for coplanar kinematics, i.e., when the ejected nucleon lies 
within the electron scattering plane P, the net ejectile polarization for an unpolarized beam 
or induced polarization (Pn), is normal to the scattering plane. In the one-photon exchange 
approximation P„ is zero when no final state interactions (FSI) between the ejected nucleon 
and the residual system are considered Thus, P„ is an observable well suited to study 

FSI effects in nuclear systems and measurements of P„ at different would give informa- 
tion about the onset of nuclear transparency. If nuclear transparency is present at certain 
value, that is, if FSI effects are quenched, we would see a decrease of P^ what would be a 
clear signature of nuclear transparency free from the ambiguities on the occupancies of the 
shells under study 0. 

The first analysis of the experiment performed at BATES by Woo et al. that measured 
Pn in ^^C(e, e'p) [jl| was made in a non relativistic framework. Other non relativistic results 
for this experiment were recently presented in reference The non relativistic approach 
to (e, e'p) is based on the impulse approximation, i.e., assuming the one-photon exchange 
picture in which the single photon interacts only with the nucleon that is detected 0. 
If no FSI are considered, the ejected nucleon is described by a plane wave (plane wave 
impulse approximation or PWIA). FSI are taken into account using potentials that distort 
the final nucleon wave function (distorted wave impulse approximation or DWIA) The 
non relativistic analyses found a a systematic underestimation of P„ of around 10% at the 
best i. 

There are two main sources of an induced nonzero normal polarization in proton knockout 
reactions caused by the interaction with the residual nucleus in the final state. One is due to 
the absorption, that is, the flux lost into inelastic channels, parametrized in the imaginary 
part of the optical potential. Semi-classically, for scattering on a given side of g, ejecting a 
nucleon from the front or the rear face of the nucleus would select out different directions of 
the angular momentum I = f x p and, as the absorption depends on how much the nucleon 
travels in the nuclear medium before being detected, the effect is a net induced polarization 
due to absorption. This is well known from hadronic reactions and is named as the Maris 
Effect or Newns Polarization |jlO|]. It is, however, a small source of P„ [||]. The bulk of 



the induced polarization is primarily due to the real part of the spin-orbit potential, that 
parametrizes the explicitly spin-dependent terms in the optical potential 

Since the spin is a property intrinsically related to relativity, one may a priori consider 
that a relativistic approach is better suited to describe nucleon polarization observables. In 
recent years, the relativistic mean-field approximation has been successfully used for the 
analyses of (e, e'p) reactions in the so-called relativistic distorted wave impulse approxima- 
tion (RDWIA) [|n|-|T^. The polarization degrees of freedom for the electron and the ejected 
nucleon have been included in this formalism years ago . In RDWIA, the nucleon current 

Jwi^, ^ = J dpippip + q)JN{^, q)i'B{p) (3) 

is calculated with relativistic ipB and ipp wave functions for initial bound and final out- 
going nucleons, respectively. is the relativistic nucleon current operator of ccl or cc2 



forms 16 . As bound state wave function, Dirac-Hartree solutions from relativistic La- 



grangian with scalar and vector (S-V) meson terms []T^ or solutions of Dirac equation with 
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phenomenological Woods-Saxon wells are customarily used. The wave function with asymp- 
totic momentum P' for the outgoing proton is a solution of the Dirac equation containing 
S-V optical potentials. Recently a relativistic calculation of P„ following those general lines 
has appeared [0 . In reference [|l8l it was found that the agreement with the data improved 
slightly compared to the non relativistic analyses of reference However, relativistic effects 
for this improvement remained unspecified in [|18] . 



Some of the differences between the relativistic and non relativistic approaches are inde- 
pendent of the dynamics, having to do with the proper (relativistic) kinematics being taken 
into account. Also, the non relativistic operators are normally obtained from an expansion 
and truncation in powers oip/M and sometimes also of g/M and uj/M. When the momenta 
and energy involved in the reaction are of the order of the nucleon mass, as it may be the 
case for (e, e'p) reactions, one must be very careful with the behavior of the expanded and 
truncated operator. In reference different non relativistic expansions were studied and 
new expressions that compared better with the unexpanded result were deduced. In refer- 
ence [pi improved non relativistic operators were used, particularly with the inclusion of the 



extra spin-orbit term in the charge density operator as described in reference [19|. This term 



proves to be necessary to reproduce at least qualitatively iT3,E0| the features seen in the 

.espce a.d TL a.„n„,et.. ... ^ ^ ^ 1" f " ^ - .neas.ed a .cent 

a[(j)' = 0) + a[(j)' = n) 

TJNAF experiment at ~ 0.8 (GeV/c)^ g. 

The non relativistic approach can be better compared to the relativistic one thinking in 
terms of the direct Pauli reduction |T^. Starting from a non relativistic formalism based 
on bispinors x solutions of a Schrodinger-like equation, one may at best construct properly 
normalized four-spinors of the form 

^-r = ^[x{p},^X{p}) (4) 

to be introduced in Eq. in order to calculate a relativistic-like nucleon current matrix 
element. In this way the relativistic kinematics is fully taken into account and no expansions 
in p/M are needed. One further step to relativize the calculations is done by rewriting the 
Dirac equation for the upper component as a Schrodinger-like equation and introducing its 
non relativistic bispinor solution x in Eq. @. This 'non relativistic' bispinor is phase-shift 
and energy eigenvalue equivalent to the relativistic solution 012^123-5^. Comparing this 



solution of the Schrodinger-like equation to the upper component of the fully relativistic 
wave functions, one finds an additional factor (i) so that the upper component of the full 
Dirac solution is quenched in the nuclear interior compared to the non relativistic solution 
1|,|2|,||,|2^. This quenching can be associated to the Darwin factor p8| that appears from 



an extra term linear in p that must be dealt with to obtain the Schrodinger-like equation 
and that is not present in the usual non relativistic treatment. 

One can then build a non relativistic formalism based on the Schrodinger-like equation, 
with central and spin-orbit potentials that are phase-shift equivalent to the relativistic po- 
tentials, incorporating a posteriori the Darwin term in order to recover exactly the same 
upper component as in RDWIA and, by means of Eq. (^), avoid the expansions in p/M. This 
formalism would incorporate all the kinematical and operator-related relativistic effects, as 
well as the dynamical quenching of the relativistic upper components due to the Darwin 
term. 
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This is done for instance by Kelly in several works p9[ though with an additional ap- 
proximation, the effective momentum approach (EMA) for the lower components. This 
amounts to approximate the a ■ p term that appears for the lower components in Eq. 

{x{p)iower = „ , ,^ X{P) upper) by ■ pas, with pas the momcutum corresponding to the 
h/ -\- ivi 

asymptotic kinematics at the nucleon vertex. Results obtained within this approximation 
both with relativistic and non relativistic potentials were compared to experiment in refer- 
ence ill. 



The differences between the calculations of [0 and those presented in |1|] can be either 
due to the EMA procedure, or to an additional dynamical relativistic effect different from the 
Darwin term, namely the enhancement of the lower components (ii): The lower components 
of the fully relativistic solutions are enhanced at the nuclear interior due to the presence of 
negative energy components [|13| , p0| , ^ . Solving the Dirac equation with scalar and vector 



potentials we see that the lower components are related to the upper ones by 

Xip)lower = ^ _|_ _|_ ^ _ yXiPjupper- (S) 

Comparing with Eq. @), we see that the lower components are enhanced with respect 

to the ones of free positive energy spinors by a factor A~^(r) = — (we 

E + M + S{r) — V{r) 

recall: 5* < 0, > and A~^{r) is ~ 2 at the nuclear interior for the usual values of the 
potentials). A~^{r) equals the inverse of the Darwin factor squared. This enhancement of 
the lower components with regard to free spinors has been sometimes referred as spinor 
distortion . 



As one can see from Eqs. (Q) and (^, for small values of the momentum p the lower 
components would play a minor role with respect to the upper ones, due to the factor 
p/ {E + M). In this low-p region {p < 300 MeV), the enhancement of the lower components 
is not important for the (e, e'p) cross sections and the most visible difference between 
RDWIA and non relativistic (kinematically corrected) DWIA results is caused by the effect 
mentioned above in point (i), namely the Darwin term. Due to dynamical effects, relativistic 
cross sections at low-p are smaller than the non relativistic ones and RDWIA-deduced 
spectroscopic factors from the \ow-p data are 10% to 15% higher than the non relativistic 
ones 



n.T^. With increasing p, however, the lower components cease to be small and their 
enhancement, present in the fully relativistic wave functions but not in Eq. (^) or in similar 
non relativistic expressions, increases the cross sections at p > 300 MeV/ c, compared to the 
non relativistic ones. This improves sizeably the agreement with the data of the RDWIA 
(e, e'p) cross sections |ll^ , |T5| . In short, in regions where the momenta of the bound and/or 
final nucleon are comparable to the nucleon mass, the RDWIA cross sections are larger 
than the non relativistic DWIA ones, in spite of the Darwin factor that in these kinematical 
regions would play a minor role. The more visible dynamical effect in high-p regions would 
be the one mentioned in paragraph (ii) that is, the enhancement of the lower components. 
This enhancement is crucial to obtain good agreement pr5| , |20[] with the recent data for Rtl 
response and Atl in ^^O taken at TJNAF for ~ 0.8 {GeV/cf |0|,|T . 



To compare with non relativistic calculations, one can project the negative energy sector 
out of the fully relativistic solutions, thus removing the enhancement of the lower compo- 
nents described in paragraph (ii). More specifically, if the negative-energy components are 
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projected out, the nucleoli current is calculated as 

Jp%(^, q) = J d0^\p + q)P{u, q)i^^^\v), (6) 

where ■?/'^-* (V'ij^'') is the positive-energy component of iI)b (V'-f), ^-c, '4'^b\p) = 
A^+){p)tljBip}, A(+)(p) = (M + J)/2M, with = {^p^ + M'\p) (similarly for ^''p^). That 
is, the matrix element of the current is computed in a truncated space with only positive 
energy spinors without enhancement of the lower components. This truncation is inherent 
to all non relativistic calculations. The dynamical enhancement of the lower components 
is contained in the current of Eq. but not in Eq. (^). Apart from kinematical effects, 
the matrix elements obtained with the prescription of Eq. are equivalent to the ones 
computed in non relativistic approaches based upon either the Foldy-Wouthuysen reduction 
@J3| or the direct Pauli reduction [0,13 . 

The EMA approach (or more properly the EMA-noSV one where, as we have said, no 
spinor distortion is considered ||2^) also removes the enhancement of the lower components 



but it is not completely equivalent to the exact projection method. Indeed, it is equivalent 
to neglecting the p-dependence of the projection operators in Eq. (^, using instead the 
asymptotic values of the momenta at the nucleon vertex. We can say that the EMA-noSV 
approach computes the matrix element with spinors that have the same structure as the ones 
that enter in the scattering of free nucleons, because it enforces the relationship between 
upper and lower components to be driven by the asymptotic value of the momenta at the 
nucleon vertex. The EMA-noSV calculation lacks any 'spinor distortion', exactly as in the 
scattering of free nucleons. In particular, the Gordon transformation is exact for the EMA- 
noSV approach and ccl and cc2 operators would lead to identical results within EMA-noSV, 
if the same choices for the off-shell values of cu, E', P and P' are made in both cases. 

The projected results, on the other hand, though lacking the large (around a factor of 
two) enhancement of the lower components seen in the fully relativistic calculation, are based 
on spinors whose upper /lower components verify Eq. (^) but with a wider value of momenta 
than in scattering from free nucleons. Thus, even projected (non relativistic) results can 
have a certain degree of spinor distortion compared to the free case due to the dispersion by 
the nuclear potentials. 

We must keep in mind that both projected and EMA-noSV results still incorporate 
the dynamical quenching of the upper components (Darwin term) and, if they are to be 
compared with non relativistic calculations, care must be taken of the Darwin term in 
the non relativistic result. Relativistic optical potentials normally give rise to increased 
absorption and stronger spin-orbit potentials. Due to this, it is expected that they would 
also lead to a stronger induced normal polarization. 

We want to emphasize that the possible differences with the former EMA-noSV analyses 
of Woo et al. |]T| are not due to the use of a relativistic optical potential. Both in references 
m and ]18[, results were presented with the same potential EDAI-C that we use in the 



present work. The Darwin term (leading to increased absorption) was also included in a 
similar way to us. Thus, if there are differences between our results and those of they 
must be due to relativistic effects additional to the Darwin term and different from the fact 
that the optical potential is relativistic or not. 
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II. RESULTS 



For the bound state wave functions we use the parameters of the set NL3 
reproduces adequately the known momentum distributions at low-p 



3311 that 



34 



We have also 

computed results with other bound state wave functions and found the effects in Pn to 
be very small up to p < 250 MeV/c. For the scattered proton wave function, we use the 
energy-dependent A-independent potentials derived by Clark et al. for ^^C (EDAI-C) and 
^^O (EDAI-0) ||35|. To study the sensitivity to different optical potentials, we also compute 
results with the energy-dependent A-dependent parametrization 1 (EDAD-1) of reference 
131. 



A. Comparison with former results from Bates 

In Fig. 1, Pn is presented against the momentum of the recoiling residual nucleus or 
missing momentum Pm, related to the momentum of the nucleon inside the target nucleus 
before being knocked-out [Q. The results are computed for the kinematics of reference [|1|, 
namely beam energy of 580 MeV, kinetic energy of the final nucleon of 270 MeV, |g| ~ 
760 MeV/c, uj ^ 290 MeV in q-uj constant kinematics with Q'^ ~ 0.5 (GeV/c)^. We use 
the Coulomb gauge in all the cases. We included the Coulomb distortion of the electron 
wave function and found its effect in P„ to be small. With solid (dotted) lines we present 
the fully relativistic results obtained with the ccl (cc2) operator. We also show results after 
projecting out the negative energy components (short-dashed lines for ccl, long-dashed lines 
for cc2). 

We see that for the fully relativistic results, the agreement with the data is excellent 
in both shells, except perhaps for the highest Pm and missing energy data points in the 
Si/2 shell, where the contribution from continuum states not considered in the present work 
begins to be important. Looking at the projected results, we see that the removal of the 
negative energy components worsens the agreement with the data for both the ccl and 
cc2 operators. In all cases P„ is smaller (less positive or more negative) for the projected 
calculations. We also see that the Gordon ambiguities, i.e., the differences between ccl 
and cc2 results, are rather small. Compared to the theoretical results of reference the 
agreement with experiment is better for any of the curves presented in Fig. 1. This cannot 
be due only to the negative energy components because the effect of projection is rather 
modest as the results in Fig. 1 show. To disentangle the reasons for this difference in Fig. 2 
we show (dotted lines) a calculation obtained with the EDAI-C potential within the EMA- 
noSV approach with the operator ccl (very similar results are obtained with the cc2 operator 
and are not shown here). As a guidance, the solid line in Fig. 2 corresponds to the same one 
of Fig. 1. We see in Fig. 2 that for the EMA-noSV results the reduction of P„ is noticeable 
and the agreement with the data is worse. Our EMA-noSV results are in the line of the ones 
obtained with the same optical potential EDAI-C in reference JH. We note that our ccl and 
cc2 results are not identical within EMA-noSV due to the different off-shell prescription we 
use in each case for the values of the kinematical quantities that enter in the evaluation of 
the current matrix element. Following reference [0 for cc2 we have 

Jtc2 = + ^^F,, (7) 
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with = {uj,q) at the electron vertex. For ccl we have 



with pm = {E,pm), Pm = P' — q, E' = y P'^ + (the final nucleon is asymptotically on- 
shell) and E = \/pm ^ + M.'^. Thus, lo used in Eq. (0) is different from uJ = E' — E a.s implied 
in Eq. (|]). If we had used U in cc2 instead of iv, the results of our ccl and cc2 EMA-noSV 
calculations would be identical. At present, there is no definite prescription for handling 
this off-shell kinematical ambiguity in u, uJ and other kinematical variables to be used in 
the current operator [23,^]. This ambiguity arises because, contrary to the scattering of 
free nucleons, part of the energy and momentum of the exchanged photon is transferred to 
the recoiling system instead of being completely absorbed by the detected nucleon. We have 
chosen the original prescription of de Forest of using E in ccl but cu in cc2. We find that, 
in this way, off-shell kinematics ambiguity effects and Gordon ambiguity ones reinforce each 
other so that the differences between the ccl and cc2 results are enhanced by our choice. 
The EMA-noSV results are free from Gordon ambiguities. This is the reason why ccl and 
cc2 results obtained in such approach are much closer than the corresponding ones of the full 
RDWIA calculations. Kinematical ambiguities will cause differences of up to 15% between 
our ccl and cc2 EMA-noSV unpolarized cross sections results for all the cases considered in 
the present work. However, the effect of these ambiguities in the P„ predictions is almost 
negligible (typically less than 1%). We have plotted only the ccl EMA-noSV result for P„. 
Our cc2 EMA-noSV result is almost identical to this curve. 

Recently, Kelly has incorporated the effect of spinor distortion within his EMA approach 
by including the relativistic potentials in the lower component of the spinors as in Eq. (^, 
while still substituting the a -p term by a-pas- This is called EMA-SV approach in reference 
29| in contrast to EMA-noSV. This procedure reintroduces the dynamical enhancement of 



the lower components. For modest values of the momentum (up to around 275 MeV/c), 
where the role of the lower components is less relevant, this approximation goes close 
to the results obtained with the exact treatment that we do in the present work. It tends, 
however, to minimize the effect of enhancement of the lower components beyond the value 
of Pm mentioned above and, as the momentum of the ejected nucleon is normally well above 
275 MeV/c, it underestimates the effect of spinor distortion for the ejected nucleon. We 
conclude that an important reason for the differences of our full RDWIA results with those 
of using the same optical potential EDAI-C, is the EMA-noSV approach employed in 
reference |1]]. 

For the purpose of comparison, we present (dash-dotted line in Fig. 2) results obtained 
with the same bound state wave functions as for the other curves, the full RDWIA approach 
with the ccl operator but a different optical potential, namely the EDAD-1 of ref. [^. We 
emphasize that the EDAI-C potential should be a more suitable choice than EDAD-1 for ^^C 
because it describes better the elastic proton scattering data for this particular nucleus. The 
effect of using EDAD-1 instead of EDAI-C is sizeable. EDAD-1 yields a larger P„ for the 
P3/2 shell, worsening the agreement with the data and a smaller P„ for the Si/2 shell, with no 
significant worse (or better) agreement with the data in this shell. We observed that both 
EDAI-C and EDAD-1 produce almost identical unpolarized cross sections or unpolarized 
response. However, the P„ values they produce are noticeably different. This shows the 
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sensitivity of Pn to details of the FSI. 

The comparison between the results of both potentials follows the same general trend as 
shown in reference [|18|. In reference [|I8|, P„ changes in the same relative direction in going 
from EDAI-C to EDAD-1 as in our calculation but, contrary to our case, the differences 
between the EDAI-C and EDAD-1 are larger for the lsi/2 shell than for the lp3/2 one and in 
this latter shell, their results seem to agree better with experiment for the EDAD-1 potential 
than for the EDAI-C. We find, in general, that our results with EDAI-C agree better with the 
data. These minor discrepancies with reference [0 should be traced back to the different 
bound state wave functions and possibly other parts of the formalism (current operator, 
Coulomb distortion) employed in both cases. 

In Fig. 2 we also show the EMA-noSV results of reference [Q with the EEI potential 
(dashed line). In reference |]1| a more positive Pn for the ps/2 shell and better agreement 
with the data was obtained with the EEI potential than with the EDAI-C one. For the lsi/2 
shell, however, a more negative P„ and worse agreement with the data was found. 

The EMA-noSV-EEI calculation of reference [|[] underestimated the data by about 10%. 
The EEI is a non relativistic optical potential obtained after folding a density-dependent 
empirical effective interaction with the nuclear density ||3^. The interaction is fitted to 
proton-nucleus elastic and inelastic scattering data for several states of several targets si- 
multaneously. The relativistic optical potentials, on the other hand, are fitted only to elastic 
scattering data. It is well known that elastic data can only constrain the asymptotic part 
of the potentials. Given the fact that the nucleus is almost transparent to electrons (com- 
pared to nucleons), phase-shift equivalent potentials that differ only in the nuclear interior 
would produce the same good fits to elastic proton scattering observables, leading however to 
different electron scattering results. The EEI approach solves this ambiguity by phenomeno- 
logically constraining the potentials in the nuclear interior by means of simultaneous fits to 
inelastic data. In the relativistic case, on the other hand, the shape of the potentials at the 
nuclear interior is assumed to be of simple Woods-Saxon + surface terms, not very differ- 
ent from what one finds in the relativistic mean field approximation. Thus, the fact that 
the relativistic model gives a very fair account of (e, e'p) observables such as P„ cannot be 
attributed merely to the incorporation of the right phenomenology, as it could be the case 
with the EEI potentials, but to a merit of the model itself. 



B. Predictions for Mainz and TJNAF in parallel kinematics 



In a recent work ||T5| , pO| , p8| , it has been shown that the dynamical enhancement of the 
lower components shows up differently in the j = 1 — 1/2 and j = 1 + 1/2 spin-orbit 
partners, specially for the Rtl response and the Atl asymmetry. We remind that for ^^C 
the two shells studied correspond to j = 1 + 1/2 spin-orbit partners. New sets of induced 
polarization P„ are being obtained at TJNAF in ^^O at a more relativistic kinematics, 
namely beam energy of 2450 MeV, kinetic energy of the ejected nucleon of about 420 MeV, 
and ~ 0.8 (GeV/c)^. There is also a proposal at Mainz to do similar measurements at 
a smaller value of the kinetic energy of the ejected nucleon, namely 200 MeV. In what follows, 
we analyze whether these experiments may provide signatures of relativistic djTiamics in P„ 
similar to the ones found in R^l and Atl- 
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In parallel kinematics only two responses, Rl and Rt, contribute to the unpolarized cross 
section and just one, R^lj Pn, so that for this kinematics the analyses are simplified. We 
present in Fig. 3 fully RDWIA results with ccl (solid line), cc2 (dotted line), projected 
results (short-dashed lines for ccl and long-dashed lines for cc2) and EMA-noSV results 
with the ccl operator (dash-dotted line) for the two p-shell spin-orbit partners of ^^O, plus 
the deep s-shell, in parallel kinematics and with beam energy and transfer energy suitable 
for Mainz ||^ (beam energy of 855 MeV, kinetic energy of the ejected nucleon of 200 MeV). 

For the three shells we see the opposite pattern to the one depicted in Fig. 1: the removal 
of the negative energy components drives here P„ towards higher values, and even more so 
does using the EMA-noSV approach. For the kinematics of Fig. 3, the effect of projection 
and Gordon ambiguities is much larger than it was in Fig. 1. 

A very characteristic feature is seen in the Si/2 shell for this case of parallel kinematics: 
A zero value of P„ is predicted within the EMA-noSV approach. A small value of P„ is 
obtained by the projected calculations, while the full RDWIA approach yields a relatively 
large (in absolute value) P„ due to spinor distortions. The choice ccl, that emphasizes the 
effect of the enhancement of the lower components |3^, yields the largest prediction for P„ in 
absolute value. Should the experiments at TJNAF or Mainz provide us with P„ values with 
equal or smaller uncertainty that the ones already measured at BATES, it will undoubtedly 
disentangle the role played in P„ by the enhancement of the lower components. 

The responses involved in the evaluation of P^ for this case are displayed in Fig 4. The 
only nonzero contribution to P„ comes from R^^^ and the link between the responses shown 
in the bottom panel of Fig. 4 and the results for P„ of Fig. 3 is straightforward. As it could be 
deduced from the values of P„ displayed in Fig. 3, R^^^ for the Si/2 shell is zero within EMA- 
noSV, it is very small for the projected results and reaches the largest absolute value for the 
full RDWIA ccl calculation. The results for the P3/2 shell follow the same trend as shown for 
the si/2 shell, only that here the more complex spin-orbit structure of the bound state causes 
a nonzero value of P„ even for the EMA-noSV results. The projected and EMA-noSV results 
display small (in absolute value) predictions for R^l- The cc2 RDWIA prediction exhibits 
larger R^l^ while the full ccl result yields the largest value of P^l- This gradation of P^^ 
is similar to what one finds generally for the unpolarized Rtl (in Q-^ constant kinematics): 
As the ccl operator enhances the role of the negative energy components ||30| with regards 
to other choices of the operator, it produces the largest value of Rtl- Thus, at least for the 
j = / -|- 1/2 spin orbit partner, we observe the same behavior for RJj^j^ and Rtl with regards 
to the effect of negative energy components. 

On the other hand, we find that projection and Gordon ambiguities effects show up 
differently for the pi/2 shell. For this case Rl, Rt and R^j^ are shown in the leftermost 
panel of Fig. 4 at the same kinematics of Fig. 3. While for RJj^j^ in the j = I + 1/2 shells the 
larger was the role given to the negative energy components the larger (in absolute value) 
RJj^L response was obtained, for the pi/2 shell (j = Z — 1/2) one sees the opposite behavior: 
the full ccl calculation yields the smallest RJ^j^ while the EMA-noSV prediction displays 
the largest one. This is at variance with the behavior observed for the unpolarized Rtl 
response and indicates an interference between positive and negative energy component 
contributions to Rtl- This interference is constructive for the j = 1 + 1/2 shells so that 
the calculations with large contribution from negative energy components yield a large Rtl^ 
while it is largely destructive for the j = / — 1/2 shells for which large effects of negative 
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energy components translate into small values of R^l- 

In reference |36| results were presented for the EMA-noSV case within the lA and also in 
a calculation beyond lA that included channel coupling to several configurations in the final 
state. Our EMA-noSV result of Fig. 4 and the one shown in Fig. 14 of reference are very 
similar, with small differences due to the different wave functions and optical potentials. The 
most interesting outcome of this comparison is that the effect of spinor distortion increases 
BJj^i, in particular for the case of the Si/2 shell that would have a zero value without spinor 
distortion within lA. Due to channel coupling (CC), a nonzero R^j^ for this shell was obtained 
in reference |PB|. The effect of spinor distortion, however, is at least twice to four times 
(depending on whether one considers the RDWIA cc2 or ccl result) larger than the one 
of CC shown in reference 



36 



We conclude that coupled channel contributions would not 
mask the large negative value of P„ caused by spinor distortion. The i?^^^. response in this 
Si/2 shell is sensitive to Gordon ambiguities and overall constitutes a very clear signature 
for the presence of negative energy components in the nucleon wave function. 

For the other shells, the effects of CC shown in reference [|3^ were small at moderate 
values of pm for the cases we studied in the present work and the lA results shown here 
should not change much if CC effects were considered. 

Still in parallel kinematics but with a larger value of that is suitable at TJNAF, we 
have obtained very similar results to the ones just presented. We plot in Fig. 5 only the 
results for P„. 



C. q-uj constant kinematics 

The experiment of Bates was performed in q-u constant kinematics and in the same 
side of q for the scattered proton, 0' = n, which corresponds to pm > in our figures. 
The analysis of this case is more complicated because all the eight responses of Eq. (0), in 
combination with the factors shown in Table 1, contribute to the cross section and Pn- In 
Fig. 6 and 7 we present the responses for the Bates results depicted in Fig. 1. The effect 
of spinor distortion and Gordon ambiguities in the Rl, Rt, Rtl and Rtt has been studied 
previously in the context of RPWIA p0| , |38| . It was found there that for the j = 1 + 1/2 



partners, as it is the case of the two shells in ^^C, the differences are relatively small, at 
least for the 'large' responses Rl and Rt- 

For the normal responses the situation is less clear. One must remember that unpolar- 
ized and normal responses share the same structure in terms of components of the hadronic 
current, differing only in the signs with which the contribution for every value of ejected 
nucleon spin projection upon the normal direction enters into the unpolarized or normal 
responses Thus, large unpolarized responses usually come from constructive interfer- 

ence of the two spin contributions and are associated with a correspondingly small normal 
response coming from destructive interference. The converse is also true: small unpolarized 
responses have a correspondingly large polarized normal response 0. If there were no FSI, 
both normal projection contributions (spin up and spin down) would be identical, all the 
responses shown in Fig. 7 would vanish and no normal polarization would be observed. 

Taking into account the value of the kinematical factors in front of each response (see 
Table I), the main contributions to -P„ for the Si/2 shell comes from RJj^ and i?^^. responses. 
In the case of the full RDWIA results, the R!},^^ response is responsible for most of the net 
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Pn- Due to this, P„ will change sign with pm because of the cos0' factor in Eq. (H). We 
also see that i?^^ is larger (in absolute value) for the calculation with larger effect of the 
negative energy components, i.e., the full ccl result. In Fig. 6 we can see the same feature in 
the unpolarized Rtl response. This characteristic has been explained before pOjpOl , ^ . For 



the p shell, the largest contribution comes from that shows little dependence on spinor 
distortion. This is why the effect of negative energy components is small for this shell. 

We plot also P„ and responses (Figs. 8, 9 and 10) for a kinematics suitable at TJNAF 
(namely |g| = 1000 MeV/c, u = 445 MeV and energy of the beam e = 2445 MeV). Apart 
from what has been already said, we find that for the P3/2 shell and Pm < there are 
small Gordon ambiguities and a very clear separation of the fully RDWIA results from the 
projected or EMA-noSV ones is seen. Therefore, this is a good region to look for the effects 
of spinor distortion. We can explain this better by looking at the results in the second 
column of Fig. 10: There, all calculations lie very close except for R^l- ^he ^3/2 shell the 
projected and EMA-noSV curves group together, while both fully relativistic calculations 
clearly deviate from the others. Going back to P„ in the second panel of Fig. 8, we observe 
these differences only in the region pm < 0, due to the different sign with which the cos (p'R^i 
term contributes in the pm < and pm > regions. This behavior is characteristic of the 
kinematics chosen at TJNAF. Indeed, as we can see in Fig. 11, at different kinematics 
conditions such as the ones suitable at Mainz, {i.e., q-uj constant kinematics with |g| = 648 
MeV/c, uj = 215 MeV and e = 855 MeV), there is not such a clear separation of the fully 
relativistic curves from the others in the P3/2 shell for pm < as the one found for the 
TJNAF kinematics. 

Another interesting feature that was already found in parallel kinematics is that, for the 
j = 1 + 1/2 shells, RJj^i has larger values when the calculation emphasizes the role of negative 
energy components while the converse is seen for the pi/2 (j = / — 1/2) shell. 

For the two p shells at the kinematics of TJNAF and Mainz, the largest contribution 
to Pn would come from VttRtt- For the P3/2 shell in the pm < region, however, this 
contribution is canceled to a large extent by the VtlRtl o^^- This explains why P„ is 
mainly negative for pm < for the P3/2 shell. On the other hand, the vtlRtl contribution 
is less important for the pi/2 shell and practically does not influence the total polarization. 
Therefore P„ for this shell is negative irrespectively of the sign of Pm- 

A serious concern is the issue of current conservation. The use of an optical potential 
breaks Gauge invariance in DWIA. We estimated the uncertainty associated to the choice 
of Gauge by comparing the results we show in the present work with the ones obtained in 
the Landau Gauge. We find that the fully relativistic results in perpendicular kinematics at 
the highest value of are the less sensitive ones to this procedure. On the other hand, the 
unpolarized cross sections can differ by as much as 50% in parallel kinematics for large values 
of \pm\- However, P„ is much less sensitive to the choice of Gauge than unpolarized cross 
sections: For the fully RDWIA results of the present work, using the Landau or Coulomb 
Gauge produces P„ results within 5%. Gauge ambiguity is much less important than the 
one due to the ccl or cc2 choice. 
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III. CONCLUSIONS 



We have found that the relativistic dynamical effect mentioned in paragraph (ii), the 
enhancement of the lower components, increases noticeably P„ with respect to both the 
projected and, more sizeably, the EMA-noSV results, driving the fully RDWIA results for 
Pn into excellent agreement with the data of reference |l|]. For the kinematics of the TJNAF 
89-033 1^ and Mainz 0] experiments we find the differences between the RDWIA and 
projected results to be important. 

Pn proves to be very sensitive to the choice of optical potential, allowing this observable 
to be used to constrain the theoretical model for FSI so that these effects can be included 
with confidence when making predictions for other observables much less sensitive to the 
choice of FSI, such as the polarization transfer observables P^ and P^ | |29| | . 

Previous explorations of the role of meson exchange currents (MEC) for Bates, based 
upon a non relativistic picture, showed very little effect in P„ at moderate Pm [§!• MEC 
are expected to play an even minor role for higher at quasi-elastic kinematics (x ~ 1) 
P] and its inclusion will not modify substantially the predictions for P„ presented in the 
present work. The same can be said of coupled channel effects analyzed within the EMA- 
noSV approach in reference |3^ . However, Gordon and kinematical off-shellness ambiguities 
are large for high experiments. We have looked for kinematical regions where these 
ambiguities are minimized. In parallel kinematics, we conclude that the Pm < region 
(|g| > \P'\) is optimal because it displays a minimum effect of Gordon ambiguities while a 
high sensitivity of the calculations to the presence of negative energy components is found. 
In q-uj constant kinematics the same favorable situation is seen again for pm < {(j)' = 0) 
but only for the P3/2 shell at |g| = 1 GeV/c, adequate to TJNAF. 

In parallel kinematics we found very clear signatures for negative energy components in 
the wave functions, that cause P„ to be driven towards more negative values with respect 
to the non relativistic prediction, particularly for the si/2 and P3/2 {j = 1 + 1/2) shells. This 
feature should remain even in the presence of MEC and CC. 

In q-u constant kinematics the effect of the negative components is manifested as an 
increase (decrease) of P„ for pm, > {pm < 0) of the relativistic predictions with regards to 
the non relativistic ones. 

Finally, we hope that future experiments will shed light on the theoretical uncertainties 
that are still present in the calculations such as which current operator should be used and 
will help to disentangle the role played by the negative energy components of the wave 
functions. 
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TABLES 



TJNAF MAMI BATES 



Oe 23.4° 48.8° 118.8° 

Q2 {GeV/cf 0.8 0.4 0.5 

oj (MeV) 445 215 292 

Vl 0.643 0.792 0.727 

Vt 0.444 0.651 3.283 

Vtl 0.737 0.932 1.642 

Vtt 0.401 0.445 0.426 



TABLE L Approximated values of the kinematical variables and the factors in Eq. (^) for the 
q-io constant experimental setups discussed in the present work. 
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FIGURES 
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FIG. 1. Pn from ^^C for tlie 1^*3/2 (upper panel) and lsi/2 (lower panel) orbits, versus missing 
momentum pm in MeV/c. Results shown correspond to a fully relativistic calculation with the 
ccl (solid) and cc2 (dotted) operators. Also shown are the results after projecting the bound and 
scattered proton wave functions over positive-energy states (short-dashed and long-dashed lines 
respectively). Data points are from j^. 
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FIG. 2. Pn from ^^C for tlie 1^*3/2 (upper panel) and lsi/2 (lower panel) orbits, versus missing 
momentum pm in MeV/c. Results shown correspond to a fully relativistic calculation with the 
ccl and the EDAI-C (solid), and the EDADl (dash-dotted) potentials. Also shown are the EMA 
results (dotted line) for the EDAI-C case. Former EMA-EEI results (dashed line) and data points 
are from 
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FIG. 3. Pn from ^^O for the lpi/2 (upper panel), 1^3/2 (mid panel) and lsi/2 (lower panel) 
orbits, versus missing momentum in MeV/c. Results shown correspond to a fully relativistic 
calculation with the ccl (solid line) and cc2 (dotted line) operators. Also shown are the projected 
results (short and long-dashed lines) and the EMA-ccl ones (dash-dotted line). Results in parallel 
kinematics corresponding to ref. and the EDAI-0 potential is used. 
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FIG. 4. Response functions corresponding to the kinematics of Fig. 3. 



19 




-300 -200 -100 100 200 300 
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FIG. 5. Pn from ^^O for the lpi/2 (upper panel), 1^3/2 (mid panel) and lsi/2 (lower panel) 
orbits, versus missing momentum in MeV/c. Results shown correspond to a fully relativistic 
calculation with the ccl (solid line) and cc2 (dotted line) operators. Also shown are the projected 
results (short and long dashed lines) and the EMA(noSV)-ccl ones (dash-dotted line). Parallel 
kinematics suitable at TJNAF and the EDAI-0 potential is used. 
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FIG. 6. Unpolarized responses for the kinematics of Fig. 1. 
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FIG. 7. Normal responses for the kinematics of Fig. 1. 
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FIG. 8. Pn from ^^O for the lpi/2 (upper panel), 1^3/2 (mid panel) and (lower panel) 

orbits, versus missing momentum pm in MeV/c. Results shown correspond to a fully relativistic 
calculation with the ccl (solid line) and cc2 (dotted line) operators. Also shown are the projected 
results (short and long dashed lines) and the EMA(noSV)-ccl ones (dash-dotted line). Results in 
q-uj constant kinematics corresponding to ref. and the EDAI-0 potential is used. 
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FIG. 9. Unpolarized responses for the kinematics of Figure 8. The response Rtt is not plotted, 
because it is much smaller than the ones shown. 
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FIG. 10. Normal responses for the kinematics of Figure 8. 
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FIG. 11. The same as Figure 8 except for slightly smaller values of q and P' as expected at 
Mainz il. 
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